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Abstract. Quantum scattering calculations are reported for the H + HC1(«, j = 0) and H + DCl(t), j — 0) 
collisions for vibrational levels v — — 2 of the diatoms. Calculations were performed for incident kinetic 
energies in the range 1CP 7 — KP 1 eV, for total angular momentum J = and s-wave scattering in 
the entrance channel of the collisions. Cross sections and rate coefficients are characterized by resonance 
structures due to quasibound states associated with the formation of the H- • -HO and H- • -DCI van der 
Waals complexes in the incident channel. For the H+HCl(v, j — 0) collision for v — 1, 2, reactive scattering 
leading to H2 formation is found to dominate over non-reactive vibrational quenching in the ultracold 
regime. Vibrational excitation of HC1 from v = to v = 2 increases the zero-temperature limiting rate 
coefficient by about 8 orders of magnitude. 

PACS. 34.50.-s Scattering of atoms and molecules - 34.50.Ez Rotational and vibrational energy transfer 
- 34. 50. Pi State-to-state scattering analyses 

1 Introduction led to the creation of molecular Bose-Einstein condensates 

starting from fermionic atoms provides unique opportu- 



Over the last several years, much progress has been made 
in cooling, trapping, and manipulating molecules at ul- 
tracold temperatures [1.23456 and Bose-Einstein con- 
densation (BEC) of diatomic molecules has recently been 
demonstrated [3|0]|3] . The experimental breakthrough that 
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nities to study the crossover regime between Bardeen- 
Cooper-Schrieffer-type superfluidity of momentum pairs 
and BEC of molecules [7H51IUIlrrj) , a topic that has been of 
long interest to the high temperature superconductivity 
community. 
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Collisional studies of ultracold molecules have received activity have not been reported. In this work, we report 

considerable attention in recent years |1 11112111. '{Ill 111 15111 61 quantum scattering calculations of atom-diatom reactions 

II 7 j and the possibility of quantum collective effects in at cold and ultracold temperatures in which the diatom 

chemical reactions involving ultracold molecules is of par- is taken to be a highly polar molecule. Previous stud- 

ticular interest |18II19| . Polar molecules are another class ies of ultracold chemical reactions investigated nonpolar 

of molecules that have received important attention in re- molecules like H2 J3U], its isotopic counterparts |4(JU41) 

cent experiments. The anisotropic, long-range character of or alkali metal dimers [I^J. Here, we investigate the bi- 

the electric dipole-dipole interactions of polar molecules molecular H + HC1 and H + DC1 collisions at low and 

also designates them as potential candidates for scalable ultralow energies for which the reaction proceeds mainly 

quantum computation schemes using electric dipole mo- by quantum tunneling of the exchanged atom through a 

ment couplings 20 21 22 23 . The techniques developed barrier along the reaction path. As elementary steps in 

so far for creating ultracold molecules fall into three dif- the H2 + CI2 reaction system, which plays a major role in 

ferent categories, namely buffer-gas cooling of paramag- chemical kinetics and in atmospheric chemistry, the gas- 

netic molecules [1.24.25,26 , electrostatic cooling of po- phase H + HC1 and CI + H2 reactions and their isotopic 

lar molecules 27,28.29.30 and photoassociation of ul- variants have received important attention both theoret- 

proaches have been successful in trapping polar molecules tally ^1.53,54,55.56.57 . However, these kinetics studies 

at cold temperatures of 10-100 mK, the creation of ultra- and experiments have been carried out in the temperature 

cold (T ~ 100 fiK) polar neutral ground state molecules range 195 K < T < 3020 K and no scattering calculations 

(KRb) was achieved only recently by photoassociation in have been reported so far for these reactions in the cold 

a magneto-optical trap |36j. Formation of electronically and ultracold regimes, to our knowledge. 

excited RbCs molecules by photoassociation in a laser- 

In this paper, we present results for state-to-state and 
cooled mixture of 85 Rb and 133 Cs atoms |3*7ll3*5) has also 

initial-state-selected probabilities and cross sections for 
been reported. 

both reactive and non-reactive channels of the H+HCl(v,j - 

Photoassociation creates molecules in highly excited 0) and H + DCl(v,j — 0) collisions for vibrational levels 

vibrational levels and their lifetime is restricted by col- v = — 2 in the ground electronic state. The presence of 

lisions leading to vibrational relaxation and/or chemical pronounced resonance structures due to quasibound states 

reactivity |39U40| . The effect of vibrational excitation on associated with van der Waals complexes in the initial 

quenching rate coefficients in the ultracold regime has channel is discussed. Finally, reaction rate coefficients for 

been explored before but similar studies on chemical re- H2 and HD formations are also presented as a function 
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of the temperature and in the zero-temperature limit. We )S 

show that vibrational excitation of HC1 and DC1 dramati- 

8e-09 

cally increases the rate coefficients in the ultracold regime. 
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2 Calculations * 4e09 

The quantum mechanical coupled-channel hypersphcrical : °""'' 
coordinate method of Skouteris et al. [HE] is used to solve 
the Schrodinger equation for the motion of the three nu- 
clei on the parametric representation of the single Born- 
Oppenheimer potential energy surface (PES) developed 
by Bian and Werner (BW) [SHj- The small effect of the 
fine-structure observed in similar reactions |52II60U61| sup- 
ports our choice to neglect the spin-orbit splitting in the 
C1( 2 P) atom. Although the accuracy of this potential en- 
ergy surface for ultracold collision studies is questionable, 
based on our experience we believe that major findings 
of our study will not be affected if a more accurate CIH2 
potential surface is used. 

Scattering calculations were performed for a total molec- 
ular angular momentum J = and s-wave scattering 
in the incident channel of the H + HC\(v,j = 0) and 
H+DCl(i>, j = 0) collisions for vibrational states v = — 2. 
We note that in the case of weak trapping potentials, 
which are expected to allow long decoherence times in 
ID trap arrays of quantum computers |23| . only s— wave 
scattering is expected to play a significant role [H21- Be- 
cause at very low kinetic and internal energies these re- 
actions proceed mainly by quantum tunneling, the result- 
ing reaction probabilities are very small and particular 
attention must be paid to convergence. Extensive con- 
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Fig. 1. Reaction probability for H2 formation in H + HCl(i> = 
0, j = 0) collisions as a function of the incident kinetic energy 
for different values of p ma x and Ap. 

vergence tests of the initial-state-selected and state-to- 
state reaction probabilities have been carried out, with 
respect to the number of rovibrational levels included in 
the basis set, j m ax, the maximum value of the hyperra- 
dius, Pmax, and the step size for the log derivative prop- 
agation, Ap. Figure n shows the total reaction probabil- 
ity for H2 formation in H + HCl(i> = 0; j = 0) collisions 
as a function of the incident kinetic energy for different 
values of p ma x and Ap. Convergence of the total reac- 
tion probability with an accuracy of the order of f CP 10 
was obtained with p max — 25.0 a.u. and Ap = 0.005 a.u. 
for kinetic energies in the range 10~ 5 — 10~ 3 eV. Using 
these values of p ma x and Ap, a more stringent conver- 
gence test consisted in the analysis of the product rota- 
tional (J') distribution. For an incident kinetic energy of 
10 -3 eV, convergence of the state-to-state probability for 
H + HCl(v = 0, 1; j = 0) and H + DCl(u = 0, f; j = 0) 
was achieved to within 10~ 10 using j max = 15 and a cut- 
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Fig. 2. State-to-state reaction probability for Hb(V = 0,j') 
formation as a function of j' in H + HC1(« = 1, j = 0) collisions, 
for different values of E max and jmax, at an incident kinetic 
energy of 10 -3 eV, pmax — 25.0 a.u., and Ap — 0.005 a.u. 
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off internal energy E max = 2.9 eV in any channel. The 
resulting basis sets for HC1 and DC1 collisions consisted of 
376 and 456 basis functions, respectively. Figure illus- 
trates the convergence of the state-to-state reaction prob- 
ability with respect to j max and E max for H2(V = Q,f) 
formation in H + HCl(w = l,j = 0) collisions, for an in- 
cident kinetic energy of 10 -3 eV, p max — 25.0 a.u., and 
Ap = 0.005 a.u. Convergence of the state-to-state proba- 
bilities for the H + HCl(t> = 2;j = 0) and H + DCl(w = 
2; j = 0) collisions was achieved using larger basis sets of 
721 and 891 basis functions, respectively, corresponding 
to jmax = 25 and E max — 3.5 eV. On the basis of these 
convergence tests, values of jmax = 15 and E max — 2.9 eV 
for v = 0, 1 and of j max = 25 and E max = 3.5 eV for v = 2 
were adopted for the calculations reported hereafter. 



Fig. 3. Cross sections for H2 and HC1 formation and nonreac- 
tive scattering in H + HCl(f, j = 0) collisions, for v — — 2, as 
a function of the incident kinetic energy. Dotted curve: nonre- 
active scattering; dashed curve: reactive HC1 product channel; 
solid curve: H2 product channel. 

3 Results and discussion 

The cross sections for H2 and HC1 formation and for non- 
reactive scattering in H + HCl(v,j = 0) collisions are 
shown in Figure |31 for v = — 2, for incident kinetic 
energy in the range 10~ 7 — lO^ 1 eV. For HC1 molecules 
initially in their ground vibrational state, the reaction pro- 
ceeds mainly by quantum tunneling through the barrier, 
yielding small values for the cross sections. We note that 
for the v = 0, j = 0, initial state nonreactive channels are 
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open only for incident kinetic energy larger than 2.585 x 
10 -3 eV, corresponding to the energy value for rotational 
excitation to the first excited level j = 1 of the ground 
vibrational state. The sharp rise in the cross section for 
nonreactive scattering at energies above 2.585 x 10 -3 eV 
is due to rotational excitation to the j = 1 level. In the 
zero-temperature limit, the cross section for the abstrac- 
tion reaction leading to H2 formation is about 9 orders of 
magnitude larger than the exchange mechanism leading 
to HC1 formation, consistent with the fact that the trans- 
mission coefficient for tunneling of the H atom through 
a finite potential barrier is larger than for heavier atoms 
like chlorine |H3]. As v increases and the barrier height 
decreases, the H2/HCI product branching ratio decreases, 
with values differing by 3 orders of magnitude for v = 1, 
and 2 orders of magnitude for v — 2 at an incident kinetic 
energy of 10~ 7 eV. For energies below I0~ 4 eV, cross sec- 
tions reach the Wigner regime |64| where they vary in- 
versely as the velocity and their ratios become constant. 
For the v = 0, j = case, the hydrogen exchange process 
is indistinguishable from elastic scattering and the cross 
section attains a constant value in the Wigner limit as 
expected for elastic scattering. For kinetic energies larger 
than 10~ 3 eV, pronounced resonance structures appear 
in the cross sections due to quasibound states associated 
with the formation of the H- • -HC1 van der Waals complex 
in the initial channel, as reported previously for different 
molecular systems {6516611671168] . Figure0|shows cross sec- 
tions for HD and HC1 formation as well as nonreactive 
scattering in H + DCl(v, j = 0) collisions as functions of 
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Fig. 4. Cross sections for HD and HC1 formation and nonre- 
active scattering in H + DCl(v,_j = 0) collisions, for v — 1, 2, as 
a function of the incident kinetic energy. Dotted curve: nonre- 
active scattering; dashed curve: reactive HC1 product channel; 
solid curve: HD product channel. 

the incident kinetic energy. Cross sections are presented 
only for the first two vibrationally excited states v = 1 
and 2 due to the negligible values obtained for the v = 
level of DC1. While nonreactive cross sections have simi- 
lar magnitude as for H + HCl(w , j = 0) collisions, reactive 
cross sections are several orders of magnitude smaller for 
the deuterated reaction, the difference being attributed to 
the less efficient tunneling of the heavier D atom. 

State-to-state reaction probabilities as a function of 
the product rotational quantum number j' are represented 
in Figure |3] for the ground vibrational state of the H2 
and HD fragments in H + HCl(u = I, j = 0) and H + 
DC1(d = I, j = 0) collisions, respectively. For a fixed in- 
cident kinetic energy of I0 -5 eV, 7 rotational levels are 
energetically accessible in the diatomic products of both 
reactions. High-j channels of H2 are preferentially popu- 
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Fig. 5. Comparison of the state-to-state reaction probabilities ^ e " 

le-16 

for H 2 (V = 0, j') (solid curve) and HD(«' = 0,/) (dot-dashed 
curve) formation in H + HCl(t> = l,j = 0) and H + DCl(w = 
l,j — 0) collisions, respectively. The probability is represented 
as a function of the product rotational number j' for a fixed 
incident kinetic energy of 10 _J eV. 
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lated even though lower-lying rotational states are open. 
The probability for H2 formation in its ground vibrational 
state peaks at j' — 4, corresponding to an exoergicity of 
0.2379 eV = 5.487 kcal/mol for the abstraction reaction. 
In the case of HD(t>' = 0) formation, the maximum energy 
released by this reaction is 0.2354 eV = 5.429 kcal/mol, 
for j' = 2, with a significantly reduced probability than 
for the H2 case. 

Figure [fj] displays the reaction rate coefficients defined 
as the product of the cross section and the relative velocity 
for H2 and HD formation for the H + HCl(t>, j = 0) and 
H + DCl(v,j = 0) reactions, respectively, for v = — 2, 
as a function of the temperature. The Wigner regime, for 
which the rate coefficients become constant, is attained 
for temperatures below 1 K for H2 formation, and below 



Fig. 6. Reaction rate coefficients for H2 and HD formation 
for the H + HCl(v, j = 0) (solid curve) and H + DC\(v,j = 0) 
(dot-dashed curve) reactions, for v = — 2, as a function of 
the temperature. 

10~ 3 K for the HD product. For v = 1 and 2, the rate 
coefficients for H2 production is an order of magnitude 
larger than for HD at cold and ultracold temperatures. 
However, for v = 1 the rate coefficient for HD formation 
becomes slightly larger than for H2 in the tunneling region 
of the H + DC1 reaction, i.e. around T = 5 K. In the zero- 
temperature limit, the rate coefficients calculated for H2 
and HD formation, respectively, are 1.9 x 10~ n cm 3 s _1 
and 1.7 x 10~ 12 cm 3 s" 1 for v = 2, 7.2 x 10~ 14 cm 3 s" 1 
and7.8xl0~ 15 cm 3 s _1 fori; = 1, and 2.4 x 10~ 19 cm 3 s" 1 
for H2 with v = 0. 
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4 Conclusion 4. M. Greiner, C. A. Regal, and D. S. Jin, Nature (London) 

426, (2003) 537. 
State-to-state and initial-state-selected cross sections have 5 M w Zwierleillj c . A . Starij c . H . Schunki s . M . Raupach) 

been calculated using fully quantum mechanical techniques S . Gupta, Z. Hafzibabic, and W. Ketterle, Phys. Rev. Lett. 

for both reactive and non-reactive channels of the H + 91, (2003) 250401. 

RC\(v, j = 0) and H + DCl(v, j = 0) collisions, for v = 6. J. Cubizolles, T. Bourdel, S. J. Kokkelmans, G. V. Shlyap- 

0-2. Resonance structures due to quasibound states of the nikov, and C. Salomon, Phys. Rev. Lett. 91, (2003) 240401. 

H- • HC1 and H- • DC1 van der Waals wells in the entrance 7 - E - Timmermans, K. Furuya, P. W. Milonni, and A. K. Ker- 

valley appear in energy dependence of the cross sections. man, Phys. Lett. A 285, (2001) 228. 

Our results also indicate that H 2 formation is the prcdom- 8 ' C ' A ' Re S aL M ' Greiner ' and D ' S ' Jin - Ph y s - Rev ' Lett ' 

inant process of H + HC1 collisions at cold and ultracold ' 



temperatures, while for H + DC1 non-reactive scattering 
is more favorable in this regime. We find that, for both 
H + HC1 and H+DC1 collisions, vibrational excitation dra- 
matically enhances the zero-temperature limiting value of 
the rate coefficients. The effect of vibrational excitation 
is found to be comparable for reactive and nonreactive 
channels in the ultracold limit. 
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